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Abstract 

The galactic population of globular clusters are old, dense star sys- 
tems, with a typical cluster containing 10 4 — 10 7 stars. As an old popu- 
lation of stars, globular clusters contain many collapsed and degenerate 
objects. As a dense population of stars, globular clusters are the scene 
of many interesting close dynamical interactions between stars. These 
dynamical interactions can alter the evolution of individual stars and 
can produce tight binary systems containing one or two compact ob- 
jects. In this review, we discuss the theoretical models of globular 
cluster evolution and binary evolution, techniques for simulating this 
evolution which lead to relativistic binaries, and current and possible 
future observational evidence for this population. Globular cluster evo- 
lution will focus on the properties that boost the production of hard 
binary systems and on the tidal interactions of the galaxy with the 
cluster, which tend to alter the structure of the globular cluster with 
time. The interaction of the components of hard binary systems alters 
the evolution of both bodies and can lead to exotic objects. Direct N- 
body integrations and Fokker-Planck simulations of the evolution of 
globular clusters that incorporate tidal interactions and lead to predic- 
tions of relativistic binary populations are also discussed. We discuss 
the current observational evidence for cataclysmic variables, millisec- 
ond pulsars, and low-mass X-ray binaries as well as possible future 
detection of relativistic binaries with gravitational radiation. 
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1 Introduction 



Relativistic binaries containing white dwarfs, neutron stars, and black holes 
in compact orbits are over-represented in globular clusters compared with 
their population in the galactic field. Observations of this population re- 
veal a host of exotic objects such as ultracompact cataclysmic variables, 
non-flickering X-ray and UV sources, low-mass X-ray binaries, and millisec- 
ond pulsars. These objects and their dark counterparts in the population 
of relativistic binaries are also likely to be observable sources of gravita- 
tional radiation for low-frequency gravitational wave detectors such as the 
planned space-borne interferometer, LISA. In the field, a relativistic binary 
is a product of the interplay between stellar evolution and the gravitational 
interaction of a tight binary. In globular clusters, the population of tight 
binaries is also a product of the dynamical evolution of an iV-body grav- 
itational system. Thus, relativistic binaries result from a combination of 
several of the more interesting processes in astrophysics. In keeping with 
the focus of this review article on relativistic binaries in globular clusters, 
we shall only touch on the aspects of globular clusters, observations, binary 
evolution, and iV-body dynamics as they relate to populations of this specific 
class of binaries in globular clusters. 

We begin by looking at the physical structure and general history of the 
galactic globular cluster system that leads to the concentration of evolved 
stars, stellar remnants, and binary systems in the cores of the clusters. 
Current observations of globular clusters that have revealed numerous pop- 
ulations of relativistic binaries and their tracers are also presented. We also 
look at the prospects for future observations in this rapidly changing area. 
Many of these relativistic binaries are the product of stellar evolution in 
compact binaries. We will look at how mass transfer from one star in the 
presence of a nearby companion can dramatically alter the evolution of both 
stars in the process of binary evolution. The enhanced production of rela- 
tivistic binaries in globular clusters results from dynamical processes that 
drive binaries toward tighter orbits and that preferentially exchange more 
massive and degenerate objects into binary systems. Numerical simulations 
of globular cluster evolution, which can be used to predict the rate at which 
relativistic binaries are formed, are discussed. These models are compared 
with the observable members of the population of relativistic binaries. Fi- 
nally, we conclude with a brief discussion of the prospects for observing these 
systems in gravitational radiation. 

Readers interested in further studies of the structure and evolution of 
globular clusters are invited to look at Binney and Tremaine |23| , Spitzer |219| , 
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and Volumes I and II of Padmanabhan's Theoretical Astrophysics 169l ll7U| 
for a good introduction to the physical processes involved. Review articles of 
Meylan & Heggie |156j and Meylan |155j also provide a comprehensive look 
at the internal dynamics of globular clusters. Although our focus is solely 
on the Galactic globular cluster system, the physics of globular cluster sys- 
tems associated with other galaxies is well covered in the review article by 
Harris [HI] as well as his lecture notes from the Saas-Fee course on star 
clusters |96j . Carney has a thorough introduction to evolution of stars in 
globular clusters |28| . An observational perspective on the role of binaries in 
globular clusters is presented in an excellent review by Bailyn while a 
good introduction to the details of observing binary systems in general can 
be found in An Introduction to Close Binary Stars jlU6j . Although slightly 
out of date, the review of binaries in globular clusters by Hut et al. |117j 
is an excellent introduction to the interaction between globular cluster dy- 
namics and binary evolution, as is a short article on globular cluster binaries 
by McMillan, Pryor, and Phinney (154] . Rappaport et al. |196j and Rasio 
et al. |197j have written reviews of numerical simulations of binary popula- 
tions in globular clusters. An excellent introduction to the astrophysics and 
numerical techniques relevant to globular cluster dynamics can be found in 
the book by Heggie and Hut [HH] . 
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2 Globular Clusters 



Globular clusters comprising 10 4 to 10 7 stars were formed early in the his- 
tory of the Milky Way, and are scattered throughout the halo. The age of 
the clusters is about 13 Gyr, with an age spread of less than 5 Gyr [23] . 
According to the frequently updated catalog of globular cluster properties 
maintained on the web by Harris |95j . the globular cluster system numbers 
roughly 150 clusters. Although the list is fairly complete, new globular clus- 
ters have been recently discovered at very low galactic latitudes |1151 1135] . 
and there is the prospect for a few more clusters to be hidden behind the 
bulge or out in the far reaches of the galaxy. The distribution of globular 
clusters in galacto-centric coordinates is shown in Figure ^ 




Figure 1: Globular cluster distribution about the galaxy. Positions are from 
Harris 195}/ and are plotted as black circles on top of the COBE FIRAS 2.2 
micron map of the Galaxy using a Mollweide projection. Figure taken from 
Brian Chaboyer's website jff^f . 

2.1 Globular cluster stars 

Because the clusters are of great age, most of the stars above about 0.8M Q 
have already evolved off the main sequence. Thus, a large number of red 
giants are readily visible in most pictures of globular clusters (see Figure |2J). 
When viewing the color-magnitude diagram (CMD) for a globular cluster, 
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Figure 2: Hubble Space Telescope photograph of the dense globular cluster 
M80 (NGC 6093). 

one can clearly see the red giant branch lifting up away from the main 
sequence. The horizontal branch of evolved stars is also seen in the CMD 
for M80 shown in Figure 

The clearly visible turn-off point in the CMD for globular clusters is 
evidence for the roughly coeval nature of the stars in the cluster. During 
the early stages of the evolution of a globular cluster, most of the gas and 
dust has been swept away. Subsequent replenishment of the intercluster gas 
by stellar winds from evolved stars is removed during periodic passages of 
the cluster through the plane of the galaxy. The remaining gas is generally 
too hot for any star formation to take place [HE]- Thus globular clusters are 
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made up of old, population II stars. 



2.2 Globular cluster structure 

The overall structure of a globular cluster can be described in terms of a 
roughly spherical iV-body system with central densities in the range 10 _1 
to 10 6 M o /pc 3 , and an average of 10 4 Mq/pc 3 . The important characteristic 
radii of a globular cluster are the core radius r c , the half-light radius r^, 
and the tidal radius r t . The core radius is defined to be the radius at which 
the surface brightness has dropped to half the central value. The half-light 
radius is the radius that contains half of the light of the cluster and the 
tidal radius is defined as the radius beyond which the external gravitational 
field of the galaxy dominates the dynamics. Theorists define to be the 
radius containing half the mass of the cluster. The half-mass radius is a 
three-dimensional theoretical construct, while the half-light radius is a two- 
dimensional observational construct. The tidal radius is always determined 
by some theoretical model. Typical values of these radii are 1.5 pc, 10 pc, 
and 50 pc, respectively |23l I170j . 

There are also important characteristic time scales that govern the dy- 
namics of globular clusters. These are the crossing time tcrossj the relaxation 
time treiaxj and the evaporation time i eva p- The crossing time is the typical 
time required for a star in the cluster to travel the characteristic size R of 
the cluster (typically taken to be the half- mass radius). Thus, t 

cross ^ 

R/v, 

where v is a typical velocity (~ 10 km/s). The relaxation time is the typical 
time for gravitational interactions with other stars in the cluster to remove 
the history of a star's original velocity. This amounts to the time required 
for gravitational encounters to alter the star's velocity by an amount com- 
parable to its original velocity. Since the relaxation time is related to the 
number and strength of the gravitational encounters of a typical cluster star, 
it is related to the number of stars in the cluster and the average energy of 
the stars in the cluster. Thus, it can be shown that the mean relaxation 
time for a cluster is [23, 169 

PAN 

^relax — ~, TT Across- (1) 

hiiV 

For a globular cluster with N = 10 5 , a characteristic size of R ~ rh ~ 10 pc, 
and a typical velocity of v ~ 10 km/s, the crossing time and relaxation 
time are t cr0 ss ~ 10 6 yr and £ re iax ~ 10 9 yr, although Binney and Tremaine 
use t cr oss ~ 10 5 yr and consequently t re i a x 10 8 yr [23] • In real globular 
clusters, the relaxation time varies throughout the cluster and the median 
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value is closer to 10 9 yr ;23_ as found in Figure 1.3 of Spitzer 219 and 
Padmanabhan |17Uj . 

The evaporation time for a cluster is the time required for the cluster 
to dissolve through the gradual loss of stars that gain sufficient velocity 
through encounters to escape its gravitational potential. In the absence of 
stellar evolution and tidal interactions with the galaxy, the evaporation time 
can be estimated by assuming that a fraction 7 of the stars in the cluster 
are evaporated every relaxation time. Thus, the rate of loss is dN/dt = 
—jN/treiax = —N/t cvap . The value of 7 can be determined by noting that 
the escape speed v e at a point x is related to the gravitational potential 
Q(x) at that point by v 2 = —2$>(x). Consequently, the mean-square escape 
speed in a cluster with density p(x) is 

2 _ Jp(x)v 2 d 3 x _ J p(x)$(x)d 3 x _ _m_ 
{ e) ~ Jp(x)d 3 x M ~ M ' { ) 

where W is the total potential energy of the cluster and M is its total mass. 
If the system is virialized (as we would expect after a relaxation time), then 
—W = 2K = M(v 2 ), where K is the total kinetic energy of the cluster, and 

(v 2 e ) = 4(v 2 ). (3) 

Thus, stars with speeds above twice the RMS speed will evaporate. Assum- 
ing a Maxwellian distribution of speeds, the fraction of stars with v > 2v Tms 
is 7 = 7.38 x 10~ 3 . Therefore, the evaporation time is 

^evap — — 136t rc i ax . (4) 

7 

Stellar evolution and tidal interactions tend to shorten the evaporation time. 
See Gnedin and Ostriker 83 and references therein for a thorough discussion 
of these effects. Using a typical t re iax for a globular cluster, we see that 
ievap ~ 10 10 yr, which is comparable to the observed age of globular clusters. 
The characteristic time scales of globular clusters differ significantly from 

each other. Across ^relax < ^ ^evap* 

When discussing stellar interactions during a given epoch of globular 
cluster evolution, it is possible to describe the background structure of the 
globular cluster in terms of a static model. These models describe the struc- 
ture of the cluster in terms of a distribution function / that can be thought 
of as providing a probability of finding a star at a particular location in 
phase-space. The static models are valid over time scales which are shorter 
than the relaxation time so that gravitational interactions do not have time 
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to significantly alter the distribution function. We can therefore assume 
df/dt ~ 0. The structure of the globular cluster is then determined by the 
collisionless Boltzmann equation, 

v • V/ - V4> ■ |£ = 0, (5) 
where the gravitational potential eft is found from / with 

V 2 </> = 4ir J f(x,v,m)d 3 vdm. (6) 

The solutions to Equation © are often described in terms of the relative 
energy per unit mass £ = — v 2 /2 with the relative potential defined as 

= —(j) + (j)Q. The constant cpQ is chosen so that there are no stars with 
relative energy less than (i.e. / > for £ > and / = for £ < 0). A 
simple class of solutions to Equation (JSJ), 

f(£) = F£ 7 / 2 , (7) 

generate what are known as Plummer models. A convenient class of models 
which admit anisotropy and a distribution in angular momenta L are known 
as King-Michie models. The King-Michie distribution function is: 



f(£, L) = Pl (2^)~ 3/2 exp [^£2) - lj , £ > 0, (8) 

with / = for £ < and p\ a constant. The velocity dispersion is de- 
termined by a and the anisotropy radius r a is defined so that the velocity 
distribution changes from nearly isotropic at the center to nearly radial at 
r a . The King-Michie distribution can be generalized to multi-mass systems, 
and although not dynamically correct, they can be used for mass estimates. 
A good description of the construction of a multi-mass King-Michie model 
can be found in the appendix of Miocchi |16L)j . 
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2.3 Globular cluster evolution 

An overview of the evolution of globular clusters can be found in Hut et 
al. |117j , Meylan & Heggie |156j , and Meylan |155j . We summarize here the 
aspects of globular cluster evolution that are relevant to the formation and 
concentration of relativistic binaries. The formation of globular clusters is 
not well understood |76| and the details of the initial mass function (IMF) 
are an ongoing field of star cluster studies. The Kroupa mass function |LS7| 



S 



is the most common IMF currently used (see |136| for a discussion of the 
local IMF). It has the form: 



dN am a% dm. (9) 

where 

a = +0.3 ± 0.7 , 0.01 < m/M Q < 0.08, 
ai = +1.3±0.5 , 0.08 < m/M e < 0.50, , , 

a 2 = +2.3 ± 0.3 , 0.50 < m/M Q < 1.00, ^ j 

a 3 = +2.3±0.7 , 1.00 < m/Af©, 

Some older work uses the Salpeter IMF which assumes a single value of a in 
Equation Q for all masses. Once the stars form out of the initial molecular 
cloud the system is not virialized (i.e. it does not satisfy Equation (jHJ), and 
it will undergo what is known as violent relaxation as the protocluster first 
begins to collapse. During violent relaxation, the total energy of individual 
stars can change as the local gravitational potential changes. The process 
of violent relaxation is a collisionless process and it occurs rapidly over a 
timescale given by a few crossing times. During violent relaxation, the 
energy per mass of a given star changes in a way that is independent of 
the mass of the star {2lT . Thus, the more massive stars will have more 
kinetic energy. These stars will then lose their kinetic energy to the less 
massive stars through stellar encounters, which leads towards equipartition 
of energy. Through virialization, this tends to concentrate the more massive 
stars in the center of the cluster — a process known as mass segregation. The 
process of mass segregation for stars of mass mj occurs on a timescale given 

by U ~ irclaxM/W 

The higher concentration of stars in the center of the cluster increases 
the probability of an encounter, which, in turn, decreases the relaxation 
time. Thus, the relaxation time given in Equation Q is an average over the 
whole cluster. The local relaxation time of the cluster is given in Meylan & 
Heggie j!56| and can be described by: 

_ 0.065(^) 3/2 

~~ p(m)G 2 In A' 1 ' 

where p is the local mass density, (v 2 ) is the mass- weighted mean square 
velocity of the stars, (m) is the mean stellar mass. The Coulomb loga- 
rithm, In A, is the logarithm of the ratio of the maximum to minimum ex- 
pected impact parameters in the cluster. Typical values of A range between 
0.4A^ [15611170] and 0.1 A?" |7S]. Binney and Tremaine provide a range of val- 
ues for In A from 10.1 in the center of the cluster to 12 at r^. Note that in the 
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central regions of the cluster, the value of t r is much lower than the average 
relaxation time. This means that in the core of the cluster, where the more 
massive stars have concentrated, there are more encounters between these 
stars. 

The concentration of massive stars in the core of the cluster will occur 
within a few relaxation times, t ~ i re iax ~ 10 9 yr. This time is longer than 
the lifetime of low metallicity stars with M > 2Mq 206 ]. Consequently, 
these stars will have evolved into carbon-oxygen (CO) and oxygen-neon 
(ONe) white dwarfs, neutron stars, and black holes. After a few more relax- 
ation times, the average mass of a star in the globular cluster will be around 
0.5M© and these degenerate objects will once again be the more massive 
objects in the cluster, despite having lost most of their mass during their 
evolution. Thus, the population in the core of the cluster will be enhanced 
in degenerate objects. Any binaries in the cluster that have a gravitational 
binding energy significantly greater than the average kinetic energy of a 
cluster star will act effectively as single objects with masses equal to their 
total mass. These objects, too, will segregate to the central regions of the 
globular cluster L 232 . The core will then be overabundant in binaries and 
degenerate objects. 

The core would undergo what is known as core collapse within a few tens 
of relaxation times unless these binaries release some of their binding energy 
to the cluster. In core collapse, the central density increases to infinity as 
the core radius shrinks to zero. An example of core collapse can be seen in 
the comparison of two cluster evolution simulations shown in Figurelll j!26j . 
Note the core collapse when the inner radius containing 1% of the total mass 
dramatically shrinks after t ~ 1 5 1 relax - Since these evolution syntheses are 
single-mass, Plummer models without binary interactions, the actual time 
of core collapse is not representative of a real globular cluster. 

The static description of the structure of globular clusters using King- 
Michie or Plummer models provides a framework for describing the envi- 
ronment of relativistic binaries and their progenitors in globular clusters. 
The short-term interactions between stars and degenerate objects can be 
analyzed in the presence of this background. Over longer time scales (com- 
parable to t re iax)) the dynamical evolution of the distribution function as well 
as population changes due to stellar evolution can alter the overall struc- 
ture of the globular cluster. We will discuss the dynamical evolution and its 
impact on relativistic binaries in Sectional 

Before moving on to the dynamical models and population syntheses of 
relativistic binaries, we will first look at the observational evidence for these 
objects in globular clusters. 
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3 Observations 



Observational evidence for relativistic binaries in globular clusters has un- 
dergone an explosion in recent years, thanks to concentrated pulsar searches, 
improved X-ray source positions from Chandra, and optical follow-ups with 
HST and ground-based telescopes. There are challenges to detecting most 
binaries since they have generally segregated to the cores of the clusters 
where crowding can be a problem. Nonetheless, numerous observations of 
both binaries and their tracer populations have been made in several glob- 
ular clusters. 

One tracer population of the dynamical processes that may lead to the 
formation of relativistic binaries is the population of blue stragglers. These 
are stars that appear on the main sequence above and to the left of the turn- 
off in the CMD of a globular cluster (see Figure EJ. These stars are hot and 
massive enough that they should have already evolved off the main sequence. 
Consequently, these objects are thought to arise from stellar coalescences 
either through the gradual merger of the components of binaries or through 
direct collisions |62l I178j . Blue stragglers are some of the most visible and 
populous evidence of the dynamical interactions that can also give rise to 
relativistic binaries. For a good description of the use of far-ultraviolet 
surveys in detecting these objects, see Knigge [132]. For somewhat older but 
still valuable reviews on the implications of blue stragglers on the dynamics 
of globular clusters, see Hut |116j and Bailyn 

The globular cluster population of white dwarfs can be used to determine 
the ages of globular clusters |162j , and so they have been the focus of targeted 
searches despite the fact that they are arguably the faintest objects in a 
globular cluster. These searches have yielded large numbers of globular 
cluster white dwarfs. For example, a recent search of u Centauri has revealed 
over 2000 white dwarfs 163 , while Hansen et al. |93| have detected 222 
white dwarfs in M4. In general, however, these searches uncover single white 
dwarfs. Optical detection of white dwarfs in binary systems tends to rely 
on properties of the accretion process related to the binary type. Therefore, 
searches for cataclysmic variables generally focus on low-luminosity X-ray 
sources |124[ 188) I229j and on ultraviolet-excess stars [Ml I133| 1151] . but 
these systems are usually a white dwarf accreting from a low mass star. The 
class of "non-flickerers" which have been detected recently |35| I224j have 
been explained as He white dwarfs in binaries containing dark CO white 
dwarfs [EEl [HH1 IH2l 

Pulsars, although easily seen in radio, are difficult to detect when they 
occur in hard binaries, due to the Doppler shift of the pulse intervals. Thanks 
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to an improved technique known as an "acceleration search" |157| . which 
assumes a constant acceleration of the pulsar during the observation period, 
more short orbital period binary pulsars are being discovered [""""J "2""] EH 
HT1 """5*1 17T1 I192| . For a good review and description of this technique, see 
Lorimer |143j . The progenitors of the ultracompact millisecond pulsars are 
thought to pass through a low-mass X-ray binary (LMXB) phase (481 1881 
I195| [198 . These systems are very bright and all of them in the globular 
cluster system are known. There are, however, several additional LMXBs 
that are currently quiescent [EB1 [230 . 

Although there are many theoretical predictions of the existence of black 
holes in globular clusters (e.g. EE3 EH US ) , there are very few ob- 
servational hints of them. Measurements of the kinematics of the cores of 
M15 [JUinn] and NGC 6752 [*"2] provide some suggestions of a large, compact 
mass. However, these observations can also be explained without requiring 
an intermediate mass black hole |148[ I173j . The unusual millisecond pulsar 
in the outskirts of NGC 6752 has also been argued to be the result of a dy- 
namical interaction with a possible binary intermediate mass black hole in 
the core |34| . If the velocity dispersion in globular clusters follows the same 
correlation to black hole mass as in galactic bulges, then there may be black 
holes with masses in the range 1 - 1O 3 M in many globular clusters |242j . 
Stellar mass black hole binaries may also be visible as low luminosity X-ray 
sources, but if they are formed in exchange interactions, they will have very 
low duty cycles and hence are unlikely to be seen |127j . 

Recent observations and catalogs of known binaries are presented in the 
following sections. 

3.1 Cataclysmic variables 

Cataclysmic variables (CVs) are white dwarfs accreting matter from a com- 
panion that is usually a dwarf star or another white dwarf. They have been 
detected in globular clusters through identification of the white dwarf it- 
self or through evidence of the accretion process. White dwarfs managed 
to avoid detection until observations with the Hubble Space Telescope re- 
vealed photometric sequences in several globular clusters j3""J ESI I172| I199| 
1201 1 12UUI 12241 19*5] . Spectral identification of white dwarfs in globular clusters 
has begun both from the ground with the VLT |161| I162j and in space with 
the Hubble Space Telescope [3*5*1 "5*51 [2"2~"fl fl63 . With spectral identification, 
it will be possible to identify those white dwarfs in hard binaries through 
Doppler shifts in the H/3 line. This approach has promise for detecting a 
large number of the expected double white dwarf binaries in globular clus- 
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ters. Photometry has also begun to reveal orbital periods |1651 1551 1128j of 
CVs in globular clusters. 

Accretion onto the white dwarf may eventually lead to a dwarf nova out- 
burst. Identifications of globular cluster CVs have been made through such 
outbursts in the cores of M5 [I5Ij . 47 Tuc Q7T], NGC 6624 |ZE|, M15 |2T2| . 
and M22 01211. With the exception of V101 in M5 j!51j . original searches 
for dwarf novae performed with ground based telescopes proved unsuccessful. 
This is primarily due to the fact that crowding obscured potential dwarf no- 
vae up to several core radii outside the center of the cluster |2()9| I211j . Since 
binaries tend to settle into the core, it is not surprising that none were found 
significantly outside of the core. Subsequent searches using the improved res- 
olution of the Hubble Space Telescope eventually revealed a few dwarf novae 
close to the cores of selected globular clusters |27)8l I2T0I 12T41 IST21 l5j . 

A more productive approach has been to look for direct evidence of the 
accretion around the white dwarf. This can be in the form of excess UV 
emission and strong Ha emission [591 1871 11331 11341 15L)| from the accretion 
disk. This technique has resulted in the discovery of candidate CVs in 
47 Tuc (2311331, M92 [EH, NGC 2808 [S3, NGC 6397 [351 EE1 E21 , and 
NGC 6712 jSU]. The accretion disk can also be discerned by very soft X- 
ray emissions. These low luminosity X-ray binaries are characterized by a 
luminosity Lx < 10 34 ' 5 erg/s, which distinguishes them from the low-mass 
X-ray binaries with Lx > 10 36 erg/s. Initial explanations of these objects 
focused on accreting white dwarfs ^U], and a significant fraction of them 
are probably CVs |230| 1233] . There have been 10 identified candidate CVs 
in 6752 H2J, 19 in 6440 QUI, 2 in u Cen [751, 5 in Terzan 5 QUI], 22 in 
47 Tuc PH, 5 in M80 |102j . and 1 in M4 [Hj. However, some of the more 
energetic sources may be LMXBs in quiescence [230j . or even candidate QSO 
sources [H] , 

The state of the field at this time is one of rapid change as Chandra 
results come in and optical counterparts are found for the new X-ray sources. 
A living catalog of CVs has been created by Downes et al. [EH an d may be 
the best source for confirmed CVs in globular clusters. 

3.2 Low- mass X-ray binaries 

The X-ray luminosities of low-mass X-ray binaries are in the range Lx ~ 
10 36 — 10 38 erg/s. The upper limit is close to the Eddington limit for ac- 
cretion onto a neutron star, so these systems must contain an accreting 
neutron star or black hole. All of the LMXBs in globular clusters contain 
an accreting neutron star as they also exhibit X-ray bursts, indicating ther- 
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monuclear flashes on the surface of the neutron star |124j . Compared with 
~ 100 such systems in the galaxy, there are 13 LMXBs known in globular 
clusters. The globular cluster system contains roughly 0.1% of the mass of 
the galaxy and roughly 10% of the LMXBs. Thus, LMXBs are substantially 
over-represented in globular clusters. 

Because these systems are so bright in X-rays, the globular cluster pop- 
ulation is completely known - we expect no new LMXBs to be discovered in 
the globular cluster system (unless more multiple sources are resolved from 
these 13 sources). The 13 sources are in 12 separate clusters. Three have 
orbital periods greater than a few hours, four ultracompact systems have 
measured orbital periods less than 1 hour, and six have undetermined orbital 
periods. The period of X1746-370 in NGC 6441 has recently been measured 
at P rb = 5.16 h using the Rossi X-ray Timing Explorer (RXTE) [12]. A 
member of the ultracompact group, 4U 1820-30 (X1820-303) in the globular 
cluster NGC 6624, has an orbital period of 11 minutes 222 . This is the 
shortest known orbital period of any binary and most certainly indicates 
a degenerate companion. The orbital period, X-ray luminosity, and host 
globular clusters for these systems are given in Tabled 



LMXB Name 


Cluster 






Ref. 






(xlO 36 erg/s) 


(hr) 




X0512-401 


NGC 1851 


1.9 


< 0.85 




X1724-307 1 


Terzan 2 


4.3 






X1730-335 


Liller 1 


2.2 




jl8ll2Tn] 


X1732-304 


Terzan 1 


0.5 




pif2T5] 


X1745-203 


NGC 6440 


0.9 






X1745-248 


Terzan 5 






EH] 


X1746-370 


NGC 6441 


7.6 


5.70 


jl8l [TTHl [2l5] 


X1747-313 


Terzan 6 


3.4 


12.36 


pi 11791 [2l5] 


X1820-303 


NGC 6624 


40.6 


0.19 


pi IT7H1 1215] 


X1832-330 


NGC 6652 


2.2 


0.73 


pilTTH] 


X1850-087 


NGC 6712 


0.8 


0.33 


pHT7mi2T5] 


X2127+119-1 


NGC 7078 


3.5 


17.10 


PI IT791 1215] 


X2127+119-2 


NGC 7078 






pi IT791 1235] 



Table 1: Low-mass X-ray binaries in globular clusters: Host clusters and 
LMXB properties. 



Sidoli et al. EE] give X1724-304. 
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The improved resolution of Chandra allows for the possibility of identify- 
ing optical counterparts to LMXBs. If an optical counterpart can be found, 
a number of additional properties and constraints for these objects can be 
determined through observations in other wavelengths. In particular, the 
orbital parameters and the nature of the secondary can be determined. So 
far, optical counterparts have been found for X0512-401 in NGC 1851 |11U| . 
X1745-203 in NGC 6440 [231], X1746-370 in NGC 6441 jUj, X1830-303 
in NGC 6624 [BT], X1832-330 in NGC 6652 [1311001, X1850-087 in NGC 
6712 [SaElIinH], X1745-248 in Terzan 5 [EI], and both LMXBs in NGC 
7078 pi I235j . Continued X-ray observations will also further elucidate the 
nature of these systems |164j . 

The 13 bright LMXBs are thought to be active members of a larger popu- 
lation of lower luminosity quiescent low mass X-ray binaries (qLMXBs) 236 . 
Early searches performed with ROSAT data (which had a detection limit 
of 10 31 erg/s) revealed roughly 30 sources in 19 globular clusters |124j . A 
more recent census of the ROSAT low luminosity X-ray sources, published 
by Verbunt |229j , lists 26 such sources that are probably related to globular 
clusters. Recent observations with the improved angular resolution of Chan- 
dra have begun to uncover numerous low luminosity X-ray candidates for 
CVs [EBJ ES2 EE3 EED EGHJ CH21 ES ESI E3 EE33 HSU • For a reasonbly com- 
plete discussion of recent observations of qLMXBs in globular clusters, see 
Verbunt and Lewin 230 or Webb and Barret |233j and references therein. 



3.3 Millisecond pulsars 

The population of known millisecond pulsars (MSPs) is one of the fastest 
growing populations of relativistic binaries in globular clusters. Several on- 
going searches are continuing to reveal millisecond pulsars in a number of 
globular clusters. Previous searches have used deep multifrequency imaging 
to estimate the population of pulsars in globular clusters (TJ. In this ap- 
proach, the expected number of pulsars beaming toward the earth, N pu \ s , 
is determined by the total radio luminosity observed when the radio beam 
width is comparable in diameter to the core of the cluster. If the minimum 
pulsar luminosity is L m [ n and the total luminosity observed is L tot , then, 
with simple assumptions on the neutron star luminosity function, 

Npuls = j 1 ^ 7T- y (12) 
-^min m { J -<tot/ - L 'min) 

In their observations of 7 globular clusters, Fruchter and Goss have recov- 
ered previously known pulsars in NGC 6440, NGC 6539, NGC 6624, and 
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47 Tuc [7^. Their estimates based on Equation (|T2*)) give evidence of a pop- 
ulation of between 60 and 200 previously unknown pulsars in Terzan 5, and 
about 15 each in Liller 1 and NGC 6544 [7T] . 

Current searches include: Arecibo, which is searching over 22 globular 
clusters j!05j : Green Bank Telescope (GBT), which is working alone and in 
conjunction with Arecibo |1231 1105] : the Giant Metrewave Radio Telescope 
(GMRT), which is searching over about 10 globular clusters ,66.; and Parkes, 
which is searching over 60 globular clusters |40j . Although these searches 
have been quite successful, they are still subject to certain selection effects 
such as distance, dispersion measure, and acceleration in compact bina- 
ries |27j . For an excellent review of the properties of all pulsars in globular 
clusters, see the review by Camilo and Rasio [2ZI and references therein. 
The properties of known pulsars in binary systems with orbital period less 
than one day are listed in Table which is a subset of Table 1 in Camilo 
and Rasio 27 a . 



Pulsar -Pspin Cluster P or b e M 2 Ref. 





(ms) 




(days) 




(M ) 




J0024-7204I 


3.484 


47 Tuc 


0.229 


< 0.0004 


0.013 




J0023-7203J 


2.100 


47 Tuc 


0.120 


< 0.00004 


0.021 


in 


J0024-7204LO 


2.643 


47 Tuc 


0.135 


< 0.00016 


0.022 


M 


J0024-72P 


3.643 


47 Tuc 


0.147 




0.017 




J0024-72R 


3.480 


47 Tuc 


0.066 




0.026 


US 


J0024-7203U 


4.342 


47 Tuc 


0.429 


0.000014 0.12 


EH 


J0024-72V 


4.810 


47 Tuc 


0.2 




0.34 


US 


J0024-7204W 


2.352 


47 Tuc 


0.133 




0.12 


P] 


J0024-72Y 


2.196 


47 Tuc 


0.521 




0.013 


EH 


B1516+02C 


2.484 


M5 


0.087 




0.037 




B1639+36D 


3.118 


M13 


0.591 




0.17 


nu 


B1639+36E 


2.487 


M13 


0.213 




0.061 




J1701-3006B 


3.593 


M62 


0.144 


< 0.00007 


0.12 


pn; 


J1701-3006C 


3.806 


M62 


0.215 


< 0.00006 


0.069 


pn; 


J1701-3006E 


3.234 


M62 


0.16 




0.03 


jSU 


J1701-3006F 


2.295 


M62 


0.20 




0.018 


jSU 


B1718-19 


1004.03 


NGC 6342 


0.258 


< 0.005 


0.11 


[22H] 


J 1748-2446 A 
J 1748- 

2446AM 
J1748-2446N 


11.563 
3.569 
8.666 


Terzan 5 
Terzan 5 
Terzan 5 


0.075 
0.443 
0.385 


< 0.0012 
0.000045 


0.087 
0.13 
0.46 


Ml 
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J 1748-24460 


1.676 


Terzan 5 


0.259 




0.035 


pun 


J1748-2446P 


1.728 


Terzan 5 


0.362 




0.36 


D322I 


J 1 748-2446 V 


2.072 


Terzan 5 


0.503 




0.11 


H22I 


J1807-24A 


3.059 


NGC 6544 


0.071 




0.0089 


HUES] 


J1911-5958A 


3.266 


NGC 6752 


0.837 


< 0.00001 


0.18 


(301 


J1911+0102A 


3.618 


NGC 6760 


0.140 


< 0.00013 


0.017 




B2127+11C 


30.529 


M 15 


0.335 


0.681 


0.92 


© 


J2140-2310A 


11.019 


M30 


0.173 


< 

0.00012 


0.01 





Table 2: Short orbital period binary millisecond pulsars in globular clusters. 
Host clusters and orbital properties. 

With the ongoing searches, it can be reasonably expected that the num- 
ber of millisecond pulsars in binary systems in globular clusters will continue 
to grow in the coming years. 



3.4 Black holes 

There have been very few observations of black hole binaries in globular 
clusters. Although there have been hints of possible black hole binaries in 
extragalactic globular cluster systems [7J > there are no known black hole 
binaries in the galactic globular cluster system. All of the globular cluster 
high luminosity LMXBs exhibit the X-ray variability that is indicative of 
nuclear burning on the surface of a neutron star. It is possible that some of 
the recently discovered low luminosity LMXBs may house black holes instead 
of neutron stars P30 , it is more likely that they are simply unusual neutron 
star LMXBs in quiesence 236 . Finally, there is very circumstantial evidence 
for the possible existence of an intermediate mass black hole (IMBH) binary 
in NGC 6752 based upon an analysis of the MSP binary PSR A [HHUHH- 
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4 Relativistic Binaries 



Relativistic binaries are binary systems with at least one degenerate or col- 
lapsed object and an orbital period such that they will be brought into 
contact within a Hubble time. (Note that this definition also includes bina- 
ries which are already in contact.) Outside of dense stellar clusters, most 
relativistic binary systems arise from primordial binary systems whose evo- 
lution drives them to tight, ultracompact orbits. The dynamical processes 
in globular clusters can drive wide binary systems toward short orbital peri- 
ods and can also insert degenerate or collapsed stars into relativistic orbits 
with other stars. Before addressing specific evolutionary scenarios, we will 
present the generic features of binary evolution that lead to the formation 
of relativistic binaries. 

4.1 Binary evolution 

The evolution of a binary system of two main-sequence (MS) stars can signif- 
icantly affect the evolution of both component stars if the orbital separation 
is sufficiently small. If the orbital period is less than about 10 days, tidal 
interactions will have circularized the orbit during the pre- and early main- 
sequence phase. [HUEnnilBO^ Both stars start in the main sequence with the 
mass of the primary M p , and the mass of the secondary M s , defined such 
that M p > M s . The binary system is described by the orbital separation 
o, and the mass ratio of the components q = M s /M p . The gravitational 
potential of the binary system is described by the Roche model where each 
star dominates the gravitational potential inside regions called Roche lobes. 
The two Roche lobes meet at the inner Lagrange point along the line joining 
the two stars. Figure El shows equipotential surfaces in the orbital plane for 
a binary with q = 0.4. If either star fills its Roche lobe, matter will stream 
from the Roche lobe filling star through the inner Lagrange point to the 
other star in a process known as Roche lobe overflow (RLOF). This mass 
transfer affects both the evolution of the components of the binary as well 
as the binary properties such as orbital period and eccentricity. 

Roche lobe overflow can be triggered by the evolution of the binary 
properties or by evolution of the component stars. On the one hand, the 
orbital separation of the binary can change so that the Roche lobe can shrink 
to within the surface of one of the stars. On the other hand, stellar evolution 
may eventually cause one of the stars to expand to fill its Roche lobe. When 
both stars in the binary are main-sequence stars, the latter process is more 
common. Since the more massive star will evolve first, it will be the first 
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to expand and fill its Roche lobe. At this stage, the mass exchange can be 
conservative (no mass is lost from the binary) or non-conservative (mass is 
lost). Depending on the details of the mass exchange and the evolutionary 
stage of the mass-losing star there are several outcomes that will lead to the 
formation of a relativistic binary. The primary star can lose its envelope, 
revealing its degenerate core as either a helium, carbon-oxygen, or oxygen- 
neon white dwarf; it can explode as a supernova, leaving behind a neutron 
star or a black hole; or it can simply lose mass to the secondary so that 
they change roles. Barring disruption of the binary, its evolution will then 
continue. In most outcomes, the secondary is now the more massive of the 
two stars and it may evolve off the main sequence to fill its Roche lobe. The 
secondary can then initiate mass transfer or mass loss with the result that 
the secondary also can become a white dwarf, neutron star, or black hole. 

The relativistic binaries that result from this process fall into a number 
of observable categories. A WD-MS or WD-WD binary may eventually be- 
come a cataclysmic variable once the white dwarf begins to accrete material 
from its companion. If the companion is a main-sequence star, RLOF can 
be triggered by the evolution of the companion. If the companion is another 
white dwarf, then RLOF is triggered by the gradual shrinking of the orbit 
through the emission of gravitational radiation. Some WD-WD cataclysmic 
variables are also known as AM CVn stars. If the total mass of the WD-WD 
binary is above the Chandrasekhar mass, the system may be a progenitor 
to a type I supernova. 

The orbit of a NS-MS or NS-WD binary will shrink due to the emission 
of gravitational radiation. At the onset of RLOF, the binary will become 
either a low-mass X-ray binary (if the donor star is a WD or MS with 
M < 2M & ), or a high- mass X-ray binary (if the donor is a more massive 
main-sequence star). These objects may further evolve to become millisec- 
ond pulsars if the NS is spun up during the X-ray binary phase |46| I198| . A 
NS-NS binary will remain virtually invisible unless one of the neutron stars 
is observable as a pulsar. A BH-MS or BH-WD binary may also become a 
low- or high-mass X-ray binary. If the neutron star is observable as a pulsar, 
a BH-NS binary will appear as a binary pulsar. BH-BH binaries will be 
invisible unless they accrete matter from the interstellar medium. A compre- 
hensive table of close binary types that can be observed in electromagnetic 
radiation can be found in Hilditch |1U6| . 

The type of binary that emerges depends upon the orbital separation 
and the masses of the component stars. During the evolution of a IOMq 
star, the radius will slowly increase by a factor of about two as the star 
progresses from zero age main sequence to terminal age main sequence. The 
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radius will then increase by about another factor of 50 as the star transitions 
to the red giant phase, and an additional factor of 10 during the transition 
to the red supergiant phase. These last two increases in size occur very 
quickly compared with the slow increase during the main-sequence evolution. 
Depending upon the orbital separation, the onset of RLOF can occur any 
time during the evolution of the star. Mass transfer can be divided into 
three cases related to the timing of the onset of RLOF. 

Case A: If the orbital separation is small enough (usually a few days), 
the star can fill its Roche lobe during its slow expansion through the 
main-sequence phase while still burning hydrogen in its core. 

Case B: If the orbital period is less than about 100 days, but longer than a 
few days, the star will fill its Roche lobe during the rapid expansion to 
a red giant with a helium core. If the helium core ignites during this 
phase and the transfer is interrupted, the mass transfer is case BB. 

Case C: If the orbital period is above 100 days, the star can evolve to the 
red supergiant phase before it fills its Roche lobe. In this case, the 
star may have a CO or ONe core. 

The typical evolution of the radius for a low metallicity star is shown in 
FigureEl Case A mass transfer occurs during the slow growth, case B during 
the first rapid expansion, and case C during the final expansion phase. The 
nature of the remnant depends upon the state of the primary during the 
onset of RLOF and the orbital properties of the resultant binary depend 
upon the details of the mass transfer. 

4.2 Mass transfer 

Although there are still many unanswered theoretical questions about the 
nature of the mass transfer phase, the basic properties of the evolution of 
a binary due to mass transfer can easily be described. The rate at which 
a star can adjust to changes in its mass is governed by three time scales. 
The dynamical time scale results from the adiabatic response of the star to 
restore hydrostatic equilibrium, and can be approximated by the free fall 
time across the radius of the star 

1/2 

min, (13) 
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where M and R are the mass and radius of the star. The thermal equilibrium 
of the star is restored over a longer period given by the thermal time scale 

GM 2 n - / M \ 2 RqLq 



RL \MqJ R L 

where L is the luminosity of the star. Finally, the main-sequence lifetime of 
the star itself provides a third time scale, which is also known as the nuclear 
time scale: 

W~7xl0 9 — -fyr. (15) 

The rate of mass transfer /loss from the Roche lobe filling star is gov- 
erned by how the star's radius changes in response to changes in its mass. 
Hjellming and Webbink jlU9j describe these changes and the response of the 
Roche lobe to mass changes in the binary using the radius-mass exponents, 
£ = din R/ din M , for each of the three processes described in Equations 1)131 
ITU IT5|) and defining 

, d In i?T 

& = < 1 + "W < 16) 

for the Roche lobe radius-mass exponent. If £l > Cdyiu the star cannot 
adjust to the Roche lobe, then the mass transfer takes place on a dynami- 
cal time scale and is limited only by the rate at which material can stream 
through the inner Lagrange point. If £dyn > Cl > Cth> then the mass trans- 
fer rate is governed by the slow expansion of the star as it relaxes toward 
thermal equilibrium, and it occurs on a thermal time scale. If both £dyn 
and Cth are greater than £l, then the mass loss is driven either by stellar 
evolution processes or by the gradual shrinkage of the orbit due to the emis- 
sion of gravitational radiation. The time scale for both of these processes is 
comparable to the nuclear time scale. A good analysis of mass transfer in 
cataclysmic variables can be found in King et al. 130 . 

Conservative mass transfer occurs when there is no mass loss from the 
system. During conservative mass transfer, the orbital elements of the binary 
can change. Consider a system with total mass M = M\ + M2 and semi- 
major axis a. The total orbital angular momentum 

1/2 

(17) 
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is a constant, and we can write a oc (M1M2) 2 . Using Kepler's third law 
and denoting the initial values by a subscript i, we find: 

P ~ [-M^\ ■ (18) 
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Differentiating Equation (|18|) and noting that conservative mass transfer 
requires M\ = —M2 gives: 



P _ 3Mi (Mi - M 2 ) 

P M X M 2 ■ [ j 

Note that if the more massive star loses mass, then the orbital period de- 
creases and the orbit shrinks. If the less massive star is the donor, then the 
orbit expands. Usually, the initial phase of RLOF takes place as the more 
massive star evolves. As a consequence, the orbit of the binary will shrink, 
driving the binary to a more compact orbit. 

In non-conservative mass transfer, both mass and angular momentum 
can be removed from the system. There are two basic non-conservative 
processes which are important to the formation of relativistic binaries - the 
common-envelope process and the supernova explosion of one component of 
the binary. The result of the first process is often a short-period, circularized 
binary containing a white dwarf. Although the most common outcome of 
the second process is the disruption of the binary, occasionally this process 
will result in an eccentric binary containing a neutron star. 

Common envelope scenarios result when one component of the binary 
expands so rapidly that the mass transfer is unstable and the companion 
becomes engulfed by the donor star. The companion then ejects the envelope 
of the donor star. The energy required to eject the envelope comes from the 
orbital energy of the binary and thus the orbit shrinks. The efficiency of 
this process determines the final orbital period after the common envelope 
phase. This is described by the efficiency parameter 

a C E = -r-= , (20) 

where AEbind is the binding energy of the mass stripped from the envelope 
and A£^ or b is the change in the orbital energy of the binary. The result of the 
process is the exposed degenerate core of the donor star in a tight, circular 
orbit with the companion. This process can result in a double degener- 
ate binary if the process is repeated twice or if the companion has already 
evolved to a white dwarf through some other process. A brief description 
of the process is outlined by Webbink |234| , and a discussion of the factors 
involved in determining cxce is presented in Sandquist et al. 205 . 

The effect on a binary of mass loss due to a supernova can be quite 
drastic. Following Padmanabhan |170j . this process is outlined using the 
example of a binary in a circular orbit with radius a. Let v be the velocity 
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of one component of the binary relative to the other component. The initial 
energy of the binary is given by 

E= l( M 1 M 2 \ 2 GM 1 M 2 _ GM1M2 m) 
2\M 1 + M 2 ) V a 2a ( ' 

Following the supernova explosion of M%, the expanding mass shell will 
quickly cross the orbit of M 2 , decreasing the gravitational force acting on 
the secondary. The new energy of the binary is then 

E , = 1 M NS M 2 2 _ GM m M 2 

where Mns is the mass of the remnant neutron star. We have assumed here 
that the passage of the mass shell by the secondary has negligible effect on 
its velocity (a safe assumption, see Pfahl et al. |175j for a discussion), and 
that the primary has received no kick from the supernova (not necessarily 
a safe assumption, but see Davies and Hansen [H2] or Pfahl et al. |176j for 
an application to globular cluster binaries). Since we have assumed that the 
instantaneous velocities of both components have not been affected, we can 
replace them by v 2 = G (Mi + M 2 ) /a, and so 

Ef= GM m M 2( M 1 + M 2 y 

2a Vm NS + M 2 J V ' 

Note that the final energy will be positive and the binary will be disrupted 
if Mns < (l/2)(Mi + M 2 ). This condition occurs when the mass ejected 
from the system is greater than half of the initial total mass: 

AM > ~ (Mi + M 2 ) , (24) 

where AM = Mi — Mns- If the binary is not disrupted, the new orbit 
becomes eccentric and expands to a new semi-major axis given by 



a' 



Mi + M 2 - AM \ 

Mi + m 2 -2Am; v ' 



and orbital period 



, ,. > 3/2 f2a'-a\ 1/2 , , 

P=P[~) —7- • (26) 



a 



We have seen that conservative mass transfer can result in a tighter bi- 
nary if the more massive star is the donor. Non-conservative mass transfer 
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can also drive the components of a binary together during a common en- 
velope phase when mass and angular momentum are lost from the system. 
Direct mass loss through a supernova explosion can also alter the properties 
of a binary, but this process generally drives the system toward larger orbital 
separation and can disrupt the binary entirely. With this exception, the im- 
portant result of all of these processes is the generation of tight binaries with 
at least one degenerate object. 

The processes discussed so far apply to the generation of relativistic bi- 
naries anywhere. They occur whenever the orbital separation of a progenitor 
binary is sufficiently small to allow for mass transfer or common envelope 
evolution. Population distributions for relativistic binaries are derived from 
an initial mass function, a distribution in mass ratios, and a distribution 
in binary separations. These initial distributions are then fed into models 
for binary evolution such as Star Track ^S] or SeBa 11891 1167j in order to 
determine rates of production of relativistic binaries. The evolution of the 
binary is often determined by the application of some simple operational 
formulae such as those described by Tout et al. |227j or Hurley et al. |112j . 
For example, Hils, Bender, and Webbink jl()8j estimated a population of 
close white dwarf binaries in the disk of the galaxy using a Salpeter mass 
function, a mass ratio distribution strongly peaked at 1, and a separation 
distribution that was flat in ln(a). Other estimates of relativistic binaries 
differ mostly by using different distributions ^3 1120| I167| I166j . 

4.3 Globular cluster processes 

When the above evolutionary scenarios are played out in the environment of 
a globular cluster, additional mechanisms arise that enhance the production 
of relativistic binaries. New binary systems can be formed by dynamical 
interaction among three or more single stars or through tidal capture, and 
the period distribution and binary components of existing binary systems 
can be altered by interactions with other stars in the cluster. We will dis- 
cuss here the broad features of these interactions and how they affect the 
evolution of binary systems toward relativistic binaries. 

The formation of binaries during the dynamical evolution of globular 
clusters can occur either through tidal capture or through iV-body interac- 
tions. Tidal capture occurs when an encounter between two stars is close 
enough that significant tides are raised on each. The tides excite non-radial 
oscillations in the stars. If the energy absorbed in these oscillations is great 
enough to leave the two stars with negative total energy, then the system 
will form a binary. This process was originally thought to be the dominant 
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channel through which binaries were formed in globular clusters |231 158j . It 
is now thought to be quite rare, as detailed calculations have shown that the 
final result is more likely to be coalescence of the two stars 11171 1198j . 
Although iV-body interactions are less likely to occur than tidally significant 
two-body interactions, they are now thought to be the dominant channel for 
the formation of binaries during the evolution of a globular cluster. This 
process, however, is not likely to produce more than a few binaries during 
the lifetime of a cluster [231 I17U| . 

Observations of present binary fractions in globular clusters combined 
with evolutionary and dynamical simulations indicate initial binary fractions 
as large as 100% are not unreasonable |121j . The existence of such a popu- 
lation of primordial binaries provides a much more efficient channel for the 
transformation of the initial distribution in component masses and orbital 
periods towards higher mass components and shorter orbital periods. This 
process follows from the interaction of primordial binaries with single stars 
and other binaries. Three results of the interaction are possible: complete 
disruption of the binary, an exchange of energy between the binary and the 
field star, or a replacement of one of the binary components by the field star. 
When a binary interacts with either a field star or with another binary, the 
energy of the interaction is shared among all stars in the interaction. The 
result is that the lowest mass object in the interaction will receive the largest 
velocity and be more likely to escape the interaction. In general, these in- 
teractions are quite complex, and must be studied numerically. A typical 
exchange interaction between a binary and a field star is shown in Figure |H1 

If the initial binding energy of the binary is large, the result of these 
interactions is to shrink the orbit of the new binary as the gravitational 
energy of the binary is used to bring the field star up to the speeds of the 
binary components. However, if the binding energy is low, the field star 
contributes energy to the components of the binary, thereby widening the 
orbit. This is an example of "Heggie's Law" [HZ], which can be summarized 
as hard binaries get harder and soft binaries get softer. For roughly equal 
mass stars, a binary is considered "hard" if its binding energy is greater 
than the average kinetic energy of a field star in the cluster and "soft" if its 
binding energy is less. For unequal mass encounters, Hills |107j has shown 
that the ratio of the orbital speeds of the binary components to the speed of 
the impactor is a better indicator of whether the binding energy will increase 
or decrease. 

The average kinetic energy of a field star in the cluster is sometimes 
related to an effective temperature of the cluster (^3 11541 1187) so that 
(mv 2 ) = 3kT. Numerical studies of the outcome of hard binary interac- 
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tions indicate that the binding energy of the binary will increase by about 
20% with each encounter |1191 1187j . Since the encounter rate is proportional 
to the semi-major axis (or 1/E) and the energy increase per encounter is 
proportional to E, the rate of hardening per relaxation time is independent 
of the energy and is (ASbind) ~ — 0.6 kT/t Te i ax (23] • A common feature of 
numerical studies of hard binary interactions is the preferential exchange of 
high-mass stars and stellar remnants with the least massive member of the 
binary |217j . Thus, the dynamical interactions in a globular cluster drive 
the initial orbital period distribution toward shorter periods by hardening 
the short period binaries while disrupting the softer binaries. Through ex- 
change interactions, the mass distribution of the binary components is also 
driven toward higher mass stars, which further enhances the number of mass- 
transferring systems that can evolve to become relativistic binaries. A very 
useful numerical simulation of multiple star interactions is Fewbody |62j . 

Because stellar remnants can also be exchanged into hard binaries, glob- 
ular cluster evolution opens up a new channel for the formation of relativistic 
binaries by introducing evolved components into binary systems that have 
not yet undergone a mass transfer phase. A particularly promising channel 
involves the exchange of a neutron star into a binary with a main-sequence 
star. The binary then undergoes case B or case C mass transfer with a 
common envelope phase, resulting in a NS-WD binary |198j . Podsiadlowski 
et al. describe a similar process without requiring the common envelope 
phase |179j . Similar interactions can occur to produce WD-WD binaries if 
a massive CO or ONe white dwarf is exchanged into a hard binary. A col- 
laboration of various groups working in stellar dynamics maintains a page 
at http://manybody.org/modest/projects.html that provides a number 



of useful computational tools for comparing how dynamical interactions can 
affect different binary evolution codes |152j . 

Black hole binaries can also form as a result of exchange interactions, but 
the process is different because black hole progenitors will evolve so quickly 
in relation to the relaxation time of most globular clusters |14()| I216j . One 
scenario that generates black hole binaries in globular clusters is described 
by Portegies Zwart and McMillan |187j . Stellar mass black holes of mass 
M ~ IOA/q will be born early in the life of a globular cluster and, through 
mass segregation, they will quickly sink to the core. Once in the core, these 
black holes will be so much more massive than the field stars that they 
will effectively form their own cluster and interact solely with themselves. 
Single black holes will form binaries with other black holes through three- 
body encounters; any black holes which are in binaries with other stars 
will team up with another black hole through exchange encounters. This 
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population of black holes and black hole binaries will then evolve separately 
from the rest of the cluster as no other stars will be massive enough to affect 
its dynamics. 

Current intermediate mass black hole (IMBH) formation scenarios that 
involve globular clusters can also affect the dynamics of the globular clus- 
ter evolution, and therefore, can affect the evolution of binaries within the 
cluster. In the two most common scenarios, an IMBH is either formed 
early in the life of the globular cluster through runaway mergers of massive 
stars |18.S( EH 0^ or it is formed through the gradual accumulation of black 
holes throughout the lifetime of the globular cluster |158j . The existence of 
an IMBH in a globular cluster can also alter its density profile, and this can 
have an affect on the rest of the dynamics of the cluster . 

We have seen how the dynamics of globular clusters can enhance the pop- 
ulation of progenitors to relativistic binaries, making the standard channels 
of mass-transfer more likely to occur. In addition, globular cluster dynamics 
can open up new channels for the formation of relativistic binaries by in- 
serting evolved, stellar remnants such as neutron stars or white dwarfs into 
binary systems and by shrinking the orbits of binary systems to enhance the 
likelihood of mass exchange. Finally, binary-single star encounters can sim- 
ply create relativistic binaries by inserting two evolved objects into a binary 
and then shrinking the orbit to ultracompact periods. We next discuss the 
probable rates for the formation of such systems and the dynamical simula- 
tions that are used to synthesize globular cluster populations of relativistic 
binaries. 
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5 Dynamical Evolution 



Simulations of the populations of relativistic binaries in globular clusters rely 
on the interplay between the evolution of individual stars in the progenitor 
systems and the evolution of globular clusters. The evolution of stars in the 
progenitor systems has been discussed in the previous section and we now 
turn to techniques for simulating the evolution of globular clusters. 

The evolution of a globular cluster is dominated by the gravitational 
interaction between the component stars in the cluster. The overall structure 
of the cluster as well as the dynamics of most of the stars in the cluster 
are determined by simple iV-body gravitational dynamics. However, the 
evolutionary time scales of stellar evolution are comparable to the relaxation 
time and core collapse time of the cluster. Consequently, stellar evolution 
affects the masses of the component stars of the cluster, which affects the 
dynamical state of the cluster. Thus, the dynamical evolution of the cluster 
is coupled to the evolutionary state of the stars. Also, as we have seen in the 
previous section, stellar evolution governs the state of the binary evolution 
and binaries may provide a means of support against core collapse. Thus, 
the details of binary evolution as coupled with stellar evolution must also be 
incorporated into any realistic model of the dynamical evolution of globular 
clusters. To close the loop, the dynamical evolution of the globular cluster 
affects the distribution and population of the binary systems in the cluster. 
In our case, we are interested in the end products of binary evolution, which 
are tied both to stellar evolution and to the dynamical evolution of the 
globular cluster. To synthesize the population of relativistic binaries, we 
need to look at the dynamical evolution of the globular cluster as well as 
the evolution of the binaries in the cluster. MODEST (MOdeling DEnse 
STellar systems), a collaboration of various groups working stellar dynamics, 
maintains a website that provides the latest information about efforts to 
combine simulations of both the dynamical evolution of iV-body systems 
and stellar evolution |152| . 

General approaches to this problem involve solving the iV-body problem 
for the component stars in the cluster and introducing binary and stellar 
evolution when appropriate to modify the iV-body evolution. There are two 
fundamental approaches to tackling this problem - direct integration of the 
equations of motion for all N bodies in the system and large- N techniques, 
such as Fokker-Planck approximations coupled with Monte Carlo treatments 
of binaries (see Heggie et al. [HEj for a comparison of these techniques). For a 
recent review of progress in implementing these techniques, see the summary 
of the MODEST-2 meeting |218| . In the next two subsections, we discuss 
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the basics of each approach and their successes and shortfalls. We conclude 
this section with a discussion of the recent relativistic binary population 
syntheses generated by dynamical simulations. 



5.1 iV-body 

The iV-body approach to modeling globular cluster dynamics involves di- 
rect calculations of the gravitational interactions between all ./V bodies in 
the simulation. In principle this is a very straightforward approach. The 
positions of the N objects in the cluster are determined by direct integration 
of the 3iV equations of motion: 

h=-y: GMj{ri ~* j) . (27) 

OH 1 J il 

When the positions indicate that objects are sufficiently close to each other, 
then the interaction is modeled to determine the outcome. In order to 
achieve realistic simulations with tidal interactions, possible mass transfer, 
and a mass spectrum of bodies, detailed stellar evolution and stellar collision 
models must be included and calculated. 

The two main codes for performing A^-body simulations are Kira and 
NBODYx. The Kira integrator is part of the Starlab environment which 
also includes stellar evolution codes and other modules for doing iV-body 
simulations |118j . The NBODYx codes have been developed and improved 
by Aarseth since the early 1960's. He has two excellent summmaries of the 
general properties and development of the NBODYx codes For further 

details, see Aarseth's book on A^-body simulations p]. A good summmary of 
general A^-body applications can also be found at the NEMO website 225 . 
NBODY6++ is a parallelization of the NBODY6 code for use on large com- 
puter clusters [220] . and a parallel version of Kira is under development 152 . 
Most large A^-body calculations are done with a special purpose computer 
called the GRAPE (GRAvity PipE) invented by Makino JTHOI - The most 
recent incarnation of the GRAPE is the GRAPE 6, which has a theoretical 
peak speed of 100 Tflops .226^ . There is also a PCI card version (GRAPE- 
6A) which is designed for use in PC clusters |72| . The GRAPE calculates the 
accelerations and jerks for the interaction between each star in the cluster. 
The next generation GRAPE-DR, which could reach ~ 1 Pflops, should be 
operational in about three years. 

The main advantage of Af-body simulations is the small number of sim- 
plifying assumptions which must be made concerning the dynamical inter- 
actions within the cluster. The specific stars and trajectories involved in 
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any interactions during the simulation are known. Therefore, the details for 
those specific interactions can be calculated during the simulation. Within 
the limits of the numerical errors that accumulate during the calculation 85 , 
one can have great confidence in the results of iV-body simulations. 

Obviously, one of the main computational difficulties is simply the CPU 
cost necessary to integrate the equations of motion for A bodies. This scales 
roughly as A 3 99 a . The other computational difficulty of the direct A-body 
method is the wide range of precision required |1161 199j . Consider the range 
of distances, from the size of neutron stars (~ 10 km) to the size of the 
globular cluster (~ 50 pc ~ 10 15 km), spanning 14 orders of magnitude. If 
the intent of the calculations is to determine the frequency of interactions 
with neutron stars, we have to know the relative position of every star to 
within 1 part in 10 14 . The range of time scales is worse yet. Considering 
that the time for a close passage of two neutron stars is on the order of 
milliseconds and that the age of a globular cluster is 10 10 yr ~ 10 20 ms, we 
find that the time scales span 20 orders of magnitude. These computational 
requirements coupled with hardware limitations mean that the number of 
bodies which can be included in a reasonable simulation is no more than 
~ 10 5 . This is about an order of magnitude less than the number of stars 
in a typical globular cluster. 

Although one has great confidence in the results of an iV-body simula- 
tion, these simulations are generally for systems that are smaller than glob- 
ular clusters. Consequently, applications of A-body simulations to globular 
cluster dynamics involve scaling lower A simulations up to the globular clus- 
ter regime. Although many processes scale with A, they do so in different 
ways. Thus, one scales the results of an A-body simulation based upon 
the assumption of a dominant process. However, one can never be certain 
that the extrapolation is smooth and that there are no critical points in the 
scaling with A. One can also scale other quantities in the model, so that 
the quantity of interest is correctly scaled |185j . An understanding of the 
nature of the scaling is crucial to understanding the applicability of A-body 
simulations to globular cluster dynamics (see Baumgardt ^3] for an exam- 
ple). The scaling problem is one of the fundamental shortcomings of the 
A-body approach. 

5.2 Fokker-Planck 

The computational limitations of A-body simulations can be sidestepped 
by describing the system in terms of distribution functions f m (x,v,t) with 
the number of stars of mass m at time t in the range (x, x + d?x) and 
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(v, v + d 3 v) given by dN = f m d 3 xd 3 v. This description requires that either 
the phase-space element d 3 xd 3 v be small enough to be infinitesimal yet 
large enough to be statistically meaningful, or that f m be interpreted as the 
probability distribution for finding a star of mass m at a location in phase 
space. The evolution of the cluster is then described by the evolution of 
f m . The gravitational interaction is provided by a smoothed gravitational 
potential <j), which is determined by 



V 2 = 4vr^ m, J / mi (x, v, t)d 3 v 



(28) 



The effect of gravitational interactions is modeled by a collision term T [f] 
(see [2*31 117U| for specific descriptions of T). The dynamics of the globular 
cluster are then governed by the Fokker-Planck equation: 

f + v-V/-Vf| = r[/]. (29) 

In the Fokker-Planck approach, the mass spectrum of stars is binned, 
with a separate f m for each bin. Increasing the resolution of the mass 
spectrum requires increasing the number of distribution functions and thus 
increasing the complexity of the problem. Consequently, Fokker-Planck 
codes can handle at most a few dozen different f m . The inclusion of addi- 
tional physical variables such as binaries adds sufficient further complexity 
that the codes are taxed beyond their capacity. Methods for numerically 
solving the Fokker-Planck equation use either an orbit-averaged form of 
Equation (H) [12], or a Monte Carlo approach p| I77 | 1781 11261 163] . 

The two time scales involved in the evolution of f m are t C ross (which 
governs changes in position) and t re \ ax (which governs changes in energy). 
The orbit-averaged form of Equation (|29fl derives from the realization that 
changes in position are essentially periodic with orbital period T ~ i C ross <C 
irelax- Thus, one can average over the rapid changes in position and retain 
the slow changes in the phase space coordinates that occur over relaxation 
times. Given suitable assumptions on the symmetry of the potential and 
the velocity distribution, when one does this the Fokker-Planck equation is 
reduced to an equation involving the energy and the magnitude of the angu- 
lar momentum. The orbit-averaged solutions of the Fokker-Planck equation 
cannot easily handle the effect of binaries and the binary interactions that oc- 
cur during the evolution of a globular cluster [73] . These effects are usually 
inserted by hand using statistical methods. The advantages of the orbit- 
averaged approach are that one can generalize it to handle anisotropy in 
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velocity, thus allowing study of the effects of the galactic gravitational field 
and tidal stripping. One can also include the rotation of the cluster |129| . 

The more recent Monte Carlo simulations |771 11261 169j do not actually 
deal with the distribution functions, but rather treat the cluster as a col- 
lection of particles that represent a spherical shell of similar stars. Based 
on the pioneering work of Henon |1U41 11U3| , they are able to represent an 
arbitrary number of species and can follow binary evolution and other ef- 
fects. The underlying treatment of relaxation throughout the simulation is 
done in the Fokker-Planck approximation, but the interactions and evolu- 
tion of the stars are handled on a particle by particle basis. Consequently, 
these codes are significantly more robust in their ability to handle realistic 
populations of stars. A nearly continuous mass spectrum can be used and 
stellar evolution and binarity can be included with relative ease. In addi- 
tion, both stellar collisions and large-angle scatterings can also be tracked. 
The primary disadvantages of these Monte Carlo codes are that they require 
spherical symmetry and that they suffer from statistical noise despite the 
large number of particles being tracked. For an excellent overview of the 
implementations and history of the Monte Carlo methods based on Henon's 
work, see Marc Freitag's link on the Working Group 3 page at the MODEST 
website |152| . 

Another approach to solving the Fokker-Planck equation makes use of 
the analogy between a globular cluster and a self-gravitating gaseous sphere |145[ 
I8U| . The most effective use of the gaseous models are in a hybrid code that 
treats the single stars in a gaseous model while treating the relaxation of 
binary, three- and four-body interactions using a Monte Carlo code |5T1I52*] . 
This approach shows promise for its flexibility in adding new physics. 

5.3 Population syntheses 

Over the last ten years, there have been several works addressing binary pop- 
ulations in globular clusters pm E2 HU El El HE E2 EH [Ell EMI EHl 
EEB21 EEE3 HSU EEHB1 EE1 EEH1 EII1 EZU . These have been derived from both 
dynamical simulations and static models. Although the motivations have 
been varied, it is often possible to extract information about the resulting 
populations of relativistic binaries. Despite the differing models and popula- 
tion synthesis techniques, the predicted populations are in rough agreement. 
Here, we summarize the different techniques and their predictions for rela- 
tivistic binaries in globular clusters. 
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5.3.1 TV-body simulations 

Although A-body simulations have the potential to provide the most de- 
tailed population syntheses of relativistic binaries in globular clusters, there 
are very few actual populations described in the literature. Most of the 
current work that treats binaries in a consistent and detailed way is lim- 
ited to open clusters |1881 11141 11131 11381 1149] and is focused on a particu- 
lar outcome of the binary population, such as blue stragglers |114j . brown 
dwarfs |138j . initial binary distributions |139| . or white dwarf CMD se- 
quences |113j . Portegies Zwart et al. focus on photometric observations of 
open clusters j!88| and on spectroscopy |184j . In their comparison of TV-body 
and Fokker-Planck simulations of the evolution of globular clusters, Taka- 
hashi and Portegies Zwart |223j followed the evolution of A = 1 K,16 K, 
and 32 K systems with initial mass functions given by Equation and 
initial density profiles set up from King models. Although they allowed for 
realistic stellar binary evolution in their comparisons, their focus was on the 
structural evolution of globular clusters. Consequently there is no binary 
population provided. Other iV-body simulations suffer from this same prob- 
lem |186| . On the other hand, recent work by Shara and Hurley has focused 
specifically on white dwarf binary populations in globular clusters and has 
produced a detailed table of close white dwarf binaries that were generated 
in their simulation |213j . 

It is possible to generate a population distribution for black hole binaries 
in globular clusters using the A-body simulations of Portegies Zwart and 
McMillan |187j that were intended to describe the population of black hole 
binaries that were ejected from globular clusters. Their scenario for black 
hole binary ejection describes a population of massive stars that evolves into 
black holes. The black holes then rapidly segregate to the core and begin 
to form binaries. As the black holes are significantly more massive than 
the other stars, they effectively form a separate sub-system, which interacts 
solely with itself. The black holes form binaries and then harden through 
binary-single black hole interactions that occasionally eject either the binary, 
the single black hole, or both. 

They simulated this scenario using A = 2048 and A = 4096 systems 
with 1% massive stars. The results of their simulations roughly confirm 
a theoretical argument based on the recoil velocity that a binary receives 
during an interaction. Noting that each encounter increases the binding 
energy by about 20% and that roughly 1/3 of this energy goes into binary 
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recoil, the minimum binding energy -Ebmin of an ejected black hole binary is 

E hmin ~ 3QW ^kT, (30) 

where (M) is the average mass of a globular cluster star and Wo = (M) \(j>o\/kT 
is the dimensionless central potential. After most binaries are ejected, 
(M) ~ 0AM & . After a few gigayears, nearly all of the black holes were 
ejected. 

At the end of this phase of black hole binary ejection, there is a 50% 
chance that a binary remains in the cluster with no other black hole to 
eject it. Thus, there should be a stellar mass black hole binary remaining in 
about half of the galactic globular clusters. The maximum binding energy of 
the remaining black hole binary is m j n and is also given by Equation 1)30(1 . 
We can then approximate the distribution in energies of the remaining black 
hole binaries as being flat in log(-Eb)- The eccentricities of this population 
will follow a thermal distribution with P(e) = 2e. 



5.3.2 Monte Carlo simulations 

Dynamical Monte Carlo simulations can be used to study the evolution of bi- 
nary populations within evolving globular cluster models. Rasio et al. |198j 
have used a Monte Carlo approach (described in Joshi et al. |125( 1126] ) to 
study the formation and evolution of NS-WD binaries, which may be pro- 
genitors of the large population of millisecond pulsars being discovered in 
globular clusters (see Section l3~3|) . In addition to producing the appropriate 
population of binary millisecond pulsars to match observations, the simu- 
lations also indicate the existence of a population of NS-WD binaries (see 
Figure EJ) . 

The tail end of the systems in group B of Figure ED represents the NS- 
WD binaries that are in very short period orbits and are undergoing a slow 
inspiral due to gravitational radiation. These few binaries can be used to 
infer an order of magnitude estimate on the population of such objects in the 
galactic globular cluster system. If we consider that there are two binaries 
with orbital period less than 2000 s out of ~ 1O 6 M0 in 47 Tuc, and assume 
that this rate is consistent throughout the globular cluster system as a whole, 
we find a total of ~ 60 such binaries. Although this estimate is quite crude, 
it compares favorably with estimates arrived at through the encounter rate 
population syntheses, which are discussed in Section [5.3.31 

More recent applications of the Monte Carlo simulations that have fo- 
cused on the properties of binaries include Fregeau et al. [H2] who look at 
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the production of blue stragglers and other collision products as a result of 
binary interactions in globular clusters and Ivanova et al. |121j who have 
studied the evolution of binary fractions in globular clusters. The latter 
work demonstrates the gradual burning of binaries in the core that delays 
the collapse of the core. In addition, they have also shown the build-up of 
short period white dwarf binaries in the core through dynamical interactions 
(see Figure ITU)) . 

There is also great promise for the hybrid gas/Monte Carlo method be- 
ing developed by Spurzem and Giersz |221| . Their recent simulation of the 
evolution of a cluster of 300,000 equal point-mass stars and 30,000 binaries 
yields a wealth of detail about the position and energy distribution of bina- 
ries in the cluster |81j . Further improvements on their code have resulted 
in direct integration of the binary-binary and binary-single interactions 82 . 
As a result, they have been able to produce empirical cross-sections for 
eccentricity variations during interactions. 

5.3.3 Encounter rate techniques 

One method for exploring the production of relativistic binary populations 
in globular clusters involves determining the encounter rate expected be- 
tween different classes of objects in a globular cluster. Sigurdsson and Phin- 
ney |217j use Monte Carlo simulations of binary encounters to infer pop- 
ulations using a static background cluster described by an isotropic King- 
Michie model. Their results are focused toward predicting the observable 
end products of binary evolution such as millisecond pulsars, cataclysmic 
variables, and blue stragglers. Therefore, there are no clear descriptions of 
relativistic binary populations provided. The work of Ivanova et al. |121j 
also uses this technique to determine the evolution of binary fractions, but 
they also do not provide sufficient detail of the population to distinguish the 
relativistic binaries from other binaries in the simulation. There is promise 
to produce a more detailed description of ultracompact X-ray binaries con- 
sisting of a white dwarf and a neutron star using encounter rates |122| . 

Davies and collaborators use the technique of calculating encounter rates 
(based on calculations of cross-sections for various binary interactions and 
number densities of stars using King-Michie static models) to determine the 
production of end products of binary evolution |45L I43j . Although they also 
do not provide a clear description of a population of relativistic binaries, 
their results allow the estimation of such a population. 

Using the encounter rates of Davies and collaborators |45[ I43j , one can 
follow the evolution of binaries injected into the core of a cluster. A frac- 
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tion of these binaries will evolve into compact binaries which will then be 
brought into contact through the emission of gravitational radiation. By 
following the evolution of these binaries from their emergence from common 
envelope to contact, we can construct a population and period distribution 
for present day globular clusters ^5] . For a globular cluster with dimension- 
less central potential Wq = 12, Davies |43| followed the evolution of 1000 
binaries over two runs. The binaries were chosen from a Salpeter IMF with 
exponent a = 2.35, and the common envelope evolution used an efficiency 
parameter uqe = 0.4. One run was terminated after 15 Gyr and the popu- 
lation of relativistic binaries which had been brought into contact through 
gravitational radiation emission was noted. The second run was allowed to 
continue until all binaries were either in merged or contact systems. There 
are four classes of relativistic binaries that are brought into contact by grav- 
itational radiation: low mass white dwarf-white dwarf binaries (WD 2 ) with 
total mass below the Chandrasekhar mass; high mass white dwarf-white 
dwarf binaries (WDj) with total mass above the Chandrasekhar mass; neu- 
tron star- white dwarf binaries (NW); and neutron star- neutron star binaries 
NS 2 . The number of systems brought into contact at the end of each run is 
given in Tabic |31 



-^evol 


WD 2 


WD 2 


NW 


NS 2 


15 Gyr 


11 





10 


1 


oo 


57 





74 


18 



Table 3: Number of relativistic binaries brought into contact through binary 
interactions. 

In the second run, the relativistic binaries had all been brought into 
contact. In similar runs, this occurs after another 15 Gyr. An estimate of 
the present-day period distribution can be made by assuming a constant 
merger rate over the second 15 Gyr. Consider the total number of binaries 
that will merge to be described by n(t). Thus, the merger rate is rj = —dn/dt. 
Assuming that the mergers are driven solely by gravitational radiation, we 
can relate n(t) to the present-day period distribution. We define n(P) to be 
the number of binaries with period less than P, and thus 

dn dn dP 

11 = ~Ht = ~dPlt' (31) 

so 

— = -=5_ f 32 ) 

dP dP/dt y ' 
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The merger rate is given by the number of mergers of each binary type 
per 1000 primordial binaries per 15 Gyr. If the orbits have been circularized 
(which is quite likely if the binaries have been formed through a common 
envelope), the evolution of the period due to gravitational radiation losses 
is given by |108| 

(W -koP~ 5/3 , (33) 



dt 



where fco is given by 



ko = ^ (2 ^)8/3iL_> [ 5/3 ! (34) 



with the "chirp mass" M 5/3 = M 1 M 2 (M 1 + M 2 )~ 1/3 . 

Following this reasoning and using the numbers in Table |3J we can deter- 
mine the present day population of relativistic binaries per 1000 primordial 
binaries. To find the population for a typical cluster, we need to determine 
the primordial binary fraction for globular clusters. Estimates of the bi- 
nary fraction in globular clusters range from 13% up to about 40% based 
on observations of either eclipsing binaries 4 , 237] I238j or luminosity func- 
tions |2031 1204] . Assuming a binary fraction of 30%, we can determine the 
number of relativistic binaries with short orbital period (P or b < -fmax) for 
a typical cluster with 1O 6 M0 and the galactic globular cluster system with 
10 Mq |217j by simply integrating the period distribution from contact P c 
up to P max 

N = [ ^-P^ 3 dP. (35) 

JPc k o 

The value of P c can be determined by using the Roche lobe radius of Eggle- 
ton [57j . 

0.49g 2 / 3 

" 0.6g 2 /3+ln(l + gi/3) a ' ^ 

and stellar radii as determined by Lynden-Bell and O'Dwyer |146j . 

The expected populations for an individual cluster and the galactic clus- 
ter system are shown in Table fusing neutron star masses of 1.4M Q , white 
dwarf masses of 0.6Mq and 0.3Mq, and -P m ax = 

2000 s. 

Although we have assumed the orbits of these binaries will be circu- 
larized, there is the possible exception of NS 2 binaries, which may have a 
thermal distribution of eccentricities if they have been formed through ex- 
change interactions rather than through a common envelope. In this case, 
Equations (|33f) and (|34|) are no longer valid. An integration over both pe- 
riod and eccentricity, using the formulae of Pierro and Pinto |177j . would be 
required. 
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Object 


WD 2 


NW 


NS 2 


Cluster 
System 


5.6 
176.5 


4.0 
125.2 


0.5 
16 



Table 4: Encounter rate estimates of the population of relativistic binaries 
in a typical globular cluster and the galactic globular cluster system. 

5.3.4 Semi-empirical methods 

The small number of observed relativistic binaries can be used to infer the 
population of dark progenitor systems For example, the low-mass X-ray 
binary systems are bright enough that we see essentially all of those that are 
in the galactic globular cluster system. If we assume that the ultracompact 
ones originate from detached WD-NS systems, then we can estimate the 
number of progenitor systems by looking at the time spent by the system 
in both phases. Let Nx be the number of ultracompact LMXBs and Tx 
be their typical lifetime. Also, let iVdet be the number of detached WD- 
NS systems that will evolve to become LMXBs, and T dct be the time spent 
during the inspiral due to the emission of gravitational radiation until the 
companion white dwarf fills its Roche lobe. If the process is stationary, we 
must have 

N x _ N dct 

-tX -tdet 

The time spent in the inspiral phase can be found from integrating Equa- 
tion (|3*3*)) to get 

^det = i~ (P 8 /3 ~ P c 8/3 ) (38) 



where Pq is the period at which the progenitor emerges from the common 
envelope and P c is the period at which RLOF from the white dwarf to 
the neutron star begins. Thus, the number of detached progenitors can be 
estimated from 

There are four known ultracompact LMXBs 48 with orbital periods 
small enough to require a degenerate white dwarf companion to the neutron 
star. There are six other LMXBs with unknown orbital periods. Thus, 4 < 
Nx < 10. The lifetime Tx is rather uncertain, depending upon the nature 
of the mass transfer and the timing when the mass transfer would cease. A 
standard treatment of mass transfer driven by gravitational radiation alone 
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gives an upper bound of Tx ~ 10 9 yr |195j . but other effects such as tidal 
heating or irradiation may shorten this to T x ~ 10 7 yr [HI HHHI - The value of 
Po depends critically upon the evolution of the neutron star-main-sequence 
binary, and is very uncertain. Both ko and P c depend upon the masses of the 
WD secondary and the NS primary. For a rough estimate, we take the mass 
of the secondary to be a typical He WD of mass OAMq and the mass of the 
primary to be 1.4M . Rather than estimate the typical period of emergence 
from the common envelope, we arbitrarily choose Po = 2000 s. We can be 
certain that all progenitors have emerged from the common envelope by the 
time the orbital period is this low. The value of P c can be determined by 
using Equation l|36jl and the radius of the white dwarf as determined by 
Lynden-Bell and O'Dwyer |146j . Adopting the optimistic values of Ax = 10 
and Tx = 10 7 yr, and evaluating Equation (|38]> gives T det ~ 10 7 yr. Thus, 
we find iVdet ~ 1 — 10, which is within an order of magnitude of the numbers 
found through dynamical simulations (Section 15.3.2(1 and encounter rate 
estimations (Section E 3.3(1 . 

Current production of ultracompact WD-NS binaries is more likely to 
arise through collisions of neutron stars with lower mass red giant stars near 
the current turn-off mass. The result of such a collision is a common enve- 
lope that will quickly eject the envelope of the red giant and leave behind 
the core in an eccentric orbit. The result of the eccentric orbit is to hasten 
the inspiral of the degenerate core into the neutron star due to gravita- 
tional radiation |174j . Consequently, T^ et can be significantly shorter 122 . 
Adopting a value of Tdet ~ 10 6 gives iVdct < 100. 

Continuing in the spirit of small number statistics, we note that there is 
one known radio pulsar in a globular cluster NS-NS binary (B2127+11C) 
and about 50 known radio pulsars in the globular cluster system as a whole 
(although this number may continue to grow) |143j . We may estimate that 
NS-NS binaries make up roughly 1/50 of the total number of neutron stars 
in the globular cluster system. A lower limit on the number of neutron stars 
comes from estimates of the total number of active radio pulsars in clusters, 
giving A r NS 2 ~ 10 5 |141j . Thus, we can estimate the total number of NS-NS 
binaries to be ~ 2000. Not all of these will be in compact orbits, but we can 
again estimate the number of systems in compact orbits by assuming that 
the systems gradually decay through gravitational radiation and thus 



where A compact is the number of systems in compact orbits (P rb < 2000 s), 
Incompact is the time spent as a compact system, and T coa i esce is the typical 



N, 



N NS 2 



compact 




(40) 
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time for a globular cluster NS-NS binary to coalesce due to gravitational 
radiation inspiral. Adopting the coalescence time of B2127+11C as typical, 
^coalesce = 2 x 10 8 yr |191j . and integrating Equation l(5B]) for two 1.4M 
neutron stars, we find -/V compac t ~ 25. Again this value compares favorably 
with the values found from encounter rate estimations. 
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6 Prospects of Gravitational Radiation 



A very exciting prospect for the observation of relativistic binaries in globu- 
lar clusters lies in the fact that they will be sources of gravitational radiation. 
There is a phase in the evolution of most relativistic binaries during which 
the orbital period is slowly shrinking due to the emission of gravitational 
radiation. If the binary is in a circularized orbit, the gravitational radiation 
will be peaked strongly in the second harmonic of the orbital period, so 
/gw = 2/ or b- Gravitational radiation can be described by the dimensionless 
strain amplitude h . Although the strength of the gravitational radiation 
varies with the orientation of the binary, an angle-averaged estimate of the 
signal strength is |142j 

At a typical globular cluster distance of r ~ 10 kpc and typical chirp mass 
of M ~ 0.5M Q , a relativistic WD-WD or WD-NS binary with P orb = 
400 s will have a gravitational wave amplitude of 10~ 22 . This value is 
within the range of the proposed space-based gravitational wave observa- 
tory lisa ing. 

Many globular clusters lie off the plane of the galaxy and are relatively 
isolated systems with known positions. The angular resolution of LISA 
improves with signal strength. By focusing the search for gravitational ra- 
diation using known positions of suspected sources, it is possible to increase 
the signal-to-noise ratio for the detected signal. Thus, the angular resolu- 
tion of LISA for globular cluster sources can be on the order of the angular 
size of the globular cluster itself at / gw > 1 mHz. Consequently, the orbital 
period distribution of a globular cluster's population of relativistic binaries 
can be determined through observations in gravitational radiation. We will 
discuss the prospects for observing each class of relativistic binaries covered 
in this review. 

White dwarf-white dwarf binaries that are formed from a common en- 
velope phase will be briefly visible while the recently revealed hot core of 
the secondary cools. These objects are most likely the "non-flickerers" of 
Cool et al. [HE] and Edmonds et al. [EE]- WD-WD binaries formed through 
exchange interactions may very well harbor white dwarfs which are too cool 
to be observed. In either case, hardening through dynamical interactions 
will become less likely as the orbit shrinks and the effective cross section of 
the binary becomes too small. These objects will then be effectively invisible 
in electromagnetic radiation until they are brought into contact and RLOF 



41 



can begin. During this invisible phase, the orbital period is ground down 
through the emission of gravitational radiation until the orbital period is a 
few hundred seconds |19|. With a frequency of 1 to 10 mHz, gravitational 
radiation from such a binary will be in the band of LISA |142j . There are 
~ 175 such systems predicted from encounter rates (see Table 0J. 

White dwarf-neutron star binaries that are expected to be progenitors 
of the millisecond pulsars must pass through a phase of gravitational radi- 
ation after the degenerate core of the donor star emerges from the common 
envelope phase and before the spin-up phase begins with the onset of mass 
transfer from the WD to the neutron star. The orbital period at the on- 
set of RLOF will be on the order of 1 to 2 minutes and the gravitational 
wave signal will be received at LISA with a signal-to-noise of 50 — 100 at 
a frequency of around 20 mHz for a globular cluster binary. Estimates of 
the number of such systems range from 1 — 10 for semi-empirical methods 
(Section l5.3.4j) to ~ 125 from encounter rates (Table 0}. 

Binaries with significant eccentricity will have a spectrum of harmonics 
of the orbital frequency, with the relative strength of the nth harmonic for 
eccentricity e given by |174j 



i 

n 



9(n,e) = - 



J n - 2 (ne) - J n _i(ne) + - J n (ne) + J n+ i(ne) - J n+2 (ne) 



n 



+(1 - e 2 ) [J n _ 2 (ne)-2J n (ne) + J n+2 (ne)] 2 + ^2 [J n (ne)] 2 |, 

(42) 

where J n is the Bessel function. The higher harmonics of sufficiently eccen- 
tric binaries (e > 0.7) can be detected by LISA even though the fundamental 
orbital frequency is well below its sensitivity band of 1 — 100 mHz |21| . 

Although the globular cluster population of NS-NS binaries is expected 
to be quite small (~ 10), they may have high eccentricities. The binary 
pulsar B2127+llCis an example of a NS-NS binary in a globular cluster. 
In terms of the unknown angle of inclination i, the companion mass to 
the pulsar is M2sini ~ 1M® and its eccentricity is e = 0.68 |143j . These 
binaries may also be detectable by LISA. If the globular cluster systems of 
other galaxies follow similar evolution as the Milky Way population, these 
binaries may be potential sources for LIGO as gravitational radiation grinds 
them down to coalescence. With their high eccentricities and large chirp 
mass, black hole binaries will also be good potential sources for gravitational 
radiation from the galactic globular cluster system |2()l I21j . 

The relatively close proximity of the galactic globular cluster system and 
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the separations between individual globular clusters allows for the identifi- 
cation of gravitational radiation sources with their individual host clusters. 
Although the expected angular resolution of LISA is not small enough to 
allow for the identification of individual sources, knowledge of the posi- 
tions of the clusters will allow for focused searches of the relativistic binary 
populations of the majority of the galactic globular clusters. Armed with 
a knowledge of the orbital periods of any detected binaries, concentrated 
searches in electromagnetic radiation can be successful in identifying rela- 
tivistic binaries that may have otherwise been missed. 
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7 Summary 



Relativistic binaries are tracers for the rich dynamical evolution of globu- 
lar clusters. The populations of these objects are the result of an interplay 
between the gravitational dynamics of large iV-body systems, the dynamics 
of mass transfer, the details of stellar evolution, and the effect of the gravi- 
tational field of the galaxy. The gravitational dynamics of globular clusters 
can enhance the population of short period binaries of main-sequence stars 
as well as inject compact objects such as white dwarfs and neutron stars 
into stellar binary systems. Once they are in such systems, the details of 
stellar evolution and mass transfer in close binary systems govern the likely 
end products of the dynamical interaction between the two stars. Further- 
more, most models of the evolution of the core of a globular cluster rely on 
the gradual hardening and ejection of binary systems to delay the onset of 
core collapse. The hardening of binaries in the core of globular clusters will 
produce relativistic binaries, but it will also eventually eject these systems 
as they gain larger and larger recoil velocities in each subsequent encounter. 
The threshold for ejection from a globular cluster depends both upon the 
gravitational potential of the cluster itself and the gravitational potential 
of its environment generated by the Milky Way. As the globular cluster 
orbits the Milky Way, its local environment changes. Consequently, if other 
dynamical processes (such as gravothermal oscillations) do not dominate, 
the globular cluster's population of relativistic binaries may also reflect the 
past orbital history of the globular cluster. 

Over the last decade, observational techniques and technology have im- 
proved to the extent that significant discoveries are being made regularly. 
At this point, the bottleneck in observations of binary millisecond pulsars, 
low-mass X-ray binaries, and cataclysmic variables is time, not technology. 
As these observational techniques are brought to bear on more clusters, 
more discoveries are bound to be made. In the next decade, the possibility 
of using gravitational wave astronomy to detect relativistic binaries brings 
the exciting possibility of identifying the populations of electromagnetically 
invisible objects such as detached WD and NS binaries and black hole bi- 
naries in globular clusters. These observations can only help to improve the 
understanding of the complex and interesting evolution of these objects and 
their host globular clusters. 
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Figure 3: Color-magnitude diagram for M80. Figure taken from the catalog 
of 52 globular clusters \202tf . The entire catalog is available at the Padova 
Globular Cluster Group website (http://dipastro.pd.astro.it/globulars/ ). 
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Figure 4: Lagrange radii indicating the evolution of a Plummer model glob- 
ular cluster for an JM -body simulation and a Monte Carlo simulation. The 
radii correspond to radii containing 0.35, 1, 3.5, 5, 7, 10, 14, 20, 30, 40, 
50, 60, 70, and 80% of the total mass. Figure taken from Joshi et al. [126]. 
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Figure 5: CMD of M3 from the Hubble Space Telescope WFPC2. Note the 
stars above and to the left of the turn-off. These are the blue stragglers. 
Figure taken from Zhao and Bailyn \24Uj - 
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Figure 6: Cross section of equipotential surfaces in the orbital plane of a 
binary with q = 0.4. The values of the potential surfaces are 5.0, 3.9075, 
3.8, 3.559, 3.2, 3.0, and 2.8. The units have been normalized to the orbital 
separation, so a = 1. 
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Figure 7: Evolution of the radius for a WMq star with Z = 0.001. Figure 
taken from Pfahl et al. \1 75\ . 
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Figure 8: Example of a binary-field star exchange interaction. The binary 
comes in from the right (red-white), while the field star (green) enters from 
the left. After a complicated interaction, the white star is ejected and the 
newly formed red-green binary is in a more tightly bound orbit. Figure taken 
from McMillan fflEJI/ . 
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Figure 9: Results of the Monte Carlo simulation of NS-WD binary gener- 
ation and evolution in ^7 Tuc. Each small dot represents a binary system. 
The circles and error bars are the 10 binary pulsars in 47 Tuc with well mea- 
sured orbits. Systems in A have evolved through mass transfer from the WD 
to the NS. Systems in B have not yet evolved through gravitational radiation 
to begin RLOF from the WD to the NS. Systems in C will not undergo a 
common envelope phase. Figure taken from Rasio et al. [1 981. 
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Figure 10: Binary period distributions from the Monte Carlo simulation of 
binary fraction evolution in ^7 Tuc. The bottom panel indicates the period 
distribution for binaries containing at least one white dwarf. N\, is the total 
number of binaries and M,, p is the number of binaries per bin. Figure taken 
from Ivanova et al. \121f . 
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